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Abstract 



Summation methods play a very important role in quantum field theory because all perturbation series are divergent 
and the expansion parameter is not always small. A number of methods have been tried in this context, most notably 
Pade approximants, Borel-Pade summation, Borel transformation with mapping, which we briefly describe and one 
on which we concentrate here, Order-Dependent Mapping (ODM). We recall the basis of the method, for a class of 
series we give intuitive arguments to explain its convergence and illustrate its properties by several simple examples. 
Since the method was proposed, some rigorous convergence proofs were given. The method has also found a number 
of applications and we shall list a few. 

Keywords: Divergent series; Summation methods; Borel transformation; Quantum mechanics. 



1. The initial motivation: Perturbative quantum field theory 

In quantum field theory, the main analytic calculation tool is the perturbative expansion. As an illustration, we 
consider the important example of the i^"* field theory [24]. In the statistical formulation, one considers the Euclidean 
(or imaginary time) action S, local functional of the field x e W^, 



where r and g are two parameters. To this action is associated a functional measure e jZ., where X is the partition 
function given by the field integral 



The limit d - Q corresponds to a simple integral. 

The case d - \ corresponds to the quantum quartic anharmonic oscillator 

Dimensions d> I correspond to quantum field theory and the expression (1) is then somewhat symbolic since the 
theory has to be modified at short distance to regularize UV divergences and renormalized to cancel them. 

In particular, the dimensions d - 2,1> are especially relevant to classical statistical physics and the theory of 
phase transitions. Finally, d = 4 is relevant to the theory of fundamental interactions at the microscopic scale. The 
corresponding relativistic quantum field theory is part of the so-called Higgs mechanism. 

For the field theory (1), the perturbative expansion amounts to an expansion in powers of the positive parameter g. 

For d> \, the difficulty of evaluating the successive perturbative terms increases very rapidly. Moreover, questions 
like regularization and renormalization arise. Therefore, the calculation of renormalization group functions in the 
d - 2> (<f>^)^ field theory up order g^ [1] is a remarkable achievement. 
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1.1. Large order behaviour of perturbative series 

In the (p'^ field theory (1), g - corresponds to a singularity since the integral (2) is not defined for g < Q. 
The perturbative series is divergent. For d < A, the large order behaviour can be inferred from a steepest descent 
calculation of the field integral (2) [19] [7]. For the quartic anharmonic oscillator {d = 1) the result was derived earlier 
from the Schrdinger equation [3]. For any physical observable /, the results have the general structure 



f, (-l)Wfc!, 



(3) 



where a depends only on d and ^ is a half-integer that depends on the observable. The coefficient A = l/a has the 
value 



d = 
d=l 
d = 2 
d = 3 



A =3/2, 

A = 8 , 

A = 35.10268957367896(1) (Zinn- Justin), 

A = 113.38350781527714(1) (Zinn- Justin). 
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(5) 
(6) 
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For d = 4,to the contribution coming from the steepest descent calculation, a contribution due to the large momentum 
singularities of Feynman diagrams has in general to be added. 

Finally, notice that for d < 4, Borel summability has been proved. 

Similar results can be obtained for a number of quantum field theories. When the formal expansion parameter 
is Planck's constant, a divergence of the form (3) is in general found (except for some fermion theories), but the 
parameter a may be complex. For an early review, see [23]. 

It follows from the large order behaviour analysis that, when the expansion parameter is not small, a summation 
of the perturbative expansion is indispensable. 



1.2. Series summation 

In the study of the fundamental interactions at the microscopic scale, it was realized that in the case of the strong 
nuclear force, unlike QED, the expansion parameter was large and, therefore, perturbation theory useless, leading 
many physicists even to reject quantum field theory as a framework to describe such phenomena. 

Before the large order behaviour was even known, in [5] it was proposed, instead, to sum the perturbative expan- 
sion, using Fade approximants and the idea was appUed to a phenomenological model, the field theory ind = 4 
dimensions. Since only two or three terms could be calculated, the possible convergence of the Fade sunmiation could 
not be checked very well. However, the results obtained in this way made much better physical sense than those of 
plain perturbation theory. For a review see [22]. 

In the seventies, one outstanding problem for which summation methods was required, is the determination of 
critical exponents and other critical quantities in the theory of second order phase transitions. Following Wilson, for 
a whole class of physical systems, these quantities can be obtained from the i4>^)^ field theory in d = 3 dimensions. 
One verifies immediately that the expansion parameter, the renormalized interaction gj, is of order 1 and a series 
summation is required (we do not discuss here the s = 4 - d expansion, but the problem is analogous). 

To deal with the practical problem of series surmnation, a method was proposed based on Borel-Pade approxi- 
mants [1]. With the knowledge of the large order behaviour, a more efficient method could be developed, combining 
a Borel transformation (actually Borel-Leroy) and a conformal mapping [17], [12], which we briefly present in next 
section. However, another method based only on the analytic properties of the series, the order-dependent mapping 
was also investigated, which we describe in more detail in section 3 (a general reference is [24]). 



2. Borel transformation and conformal mapping 

The values of critical exponents in a large class of continuous (or second order) phase transitions can be inferred 
from so-called renormalization group (RG) functions of the (0^)^ quantum field theory. One important function is 
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Table 1: 

Series summed by the method based on Borel transformation and mapping for the zero g* of the RG P(g) function and 

the exponents y and v in the <p'^ field theory. 



k 


2 


3 


4 


5 


6 


7 


8* 


1.8774 


1.5135 


1.4149 


1.4107 


1.4103 


1.4105 


V 


0.6338 


0.6328 


0.62966 


0.6302 


0.6302 


0.6302 


y 


1.2257 


1.2370 


1.2386 


1.2398 


1.2398 


1.2398 



the RG y6-function whose zeros determine the RG fixed points. For exiimple, in the case of the 0* theory in J = 3 
dimensions, Nickel [1] has calculated 

308 

^(1) = -I + 0.35106959771'' -0.37652682831^ +0.495547511^ 

-0.7496891^ + (8) 

where g - 3gr/(16;r) and gi is the so-called renormalized interaction, related to the parameter that appears in the 
action (1) by = ^ + Oig^) and 

Agr) = -^AI) = - ,, • (9) 

3 dln^/d^r 

The perturbative expansion is divergent (equation (3)). For the jS-function in three dimensions, 

k 

the large order behaviour, implied by the estimate (7), is given by 

j8i oc {-afk'l^kl 

k—*co 

withfl = 0.147774232.... 

To characterize the large distance properties of statistical systems at the phase transition, one must first determine 
the non-trivial zero g* of the yS-function and then calculate various physical quantities like critical exponents f or | = |* . 
One discovers that g* is a number of order 1 and, thus, a numerical determination from the series 8 clearly requires a 
summation of the series. 

In three dimensions, the perturbative expansion is proved to be Borel summable. It is thus natural to introduce the 
Borel-Laplace transformation (here, Borel-Leroy): 

where <t is a free parameter. Then, formally in the sense of power series 



The function B^ig) is analytic in a circle of radius 1/a. The series is said Borel summable if, in addition, Bg-ig) is 
analytic in a neighbourhood of the real positive semi-axis and the integral converges. 

The series defines the function in a circle. It is thus necessary to perform an analytic continuation. In practice, 
with a small number of terms, the continuation requires a domain of analyticity larger than rigorously established. 



Le Guillou and Ziiin- Justin [17] have assumed maximal analyticity, i.e., analyticity in a cut-plane. The continuation 
has then be obtained by a conformal mapping of the cut-plane onto a circle. Finally, various modifications has been 
introduced to optimize the summation method (for details see [17]). 

Further optimization of the summation technique and the additional seven-loop contributions have led to new 
estimates of critical exponents [12]. Some results are displayed in table 1. 

3. Order-dependent mapping 

The order-dependent mapping (ODM) summation method [20] is based on some knowledge of the analytic prop- 
erties of the fimction that is expanded. It applies both to convergent and divergent series, although it is mainly useful 
in the latter case. 

3.1. The general method 

Let f{z) be an analytic function that has the Taylor series expansion 

(the = sign has to understood in the sense of series expansion.) 

When the Taylor series has a finite radius of convergence, to continue the function in the whole domain of analyt- 
icity, one can map the domain onto a circle, while preserving the origin. 

Divergent series: the intuitive idea. In a case of a divergent series, one adds to the domain of analyticity a disk 
\z\ < r of variable radius r and applies a similar mapping. Of course, the transformed series is still divergent. Then, 
one recalls the empirical rule that, for a divergent series, one is instructed to truncate the series at the term of minimal 
modulus, the last term giving an order of magnitude of the error. By adjusting the radius r order by order, one can 
manage to set the minimum always just at the last calculated term. 
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Figure 1: Mapping z i-» /I: example of a function analytic in a cut-plane. 

In what follows, we consider only functions analytic in a sector (as in the example of figure 1) and mappings 
z i-» A of the form 

z = pf(A ^(i) = A + <9(i2), 

where ^(/l) is an explicit analytic function and p an adjustable parameter. 

Although the transformed series is still divergent at p fixed, we shall verify on a few examples that, by adjusting p 
order by order (here, we hmit ourselves to Borel summable examples) one can construct a convergent algorithm. 

After the transformation, / is given by a Taylor series in A of the form 

fiz(A)) = J]Pk(p)A', 

k=0 
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where the coefficients Pkip) are polynomials of degree k in p. Since the result is formally independent of the parameter 
p, the parameter can be chosen freely. 

The ^-th approximant /**'(z) is constructed in the following way: one truncates the expansion at order k and 
chooses p as to cancel the last term. Since Pk(p) has k roots (real or complex), one chooses for p the largest possible 
root (in modulus) pk for which P'j^ip) is small. This leads to a sequence of approximants 



In the case of convergent series, it is expected that pk has a non-vanishing limit for A: — > oo. By contrast, for divergent 
series it is expected that p^ goes to zero for large k as 



The intuitive idea here is that pk corresponds to a 'local' radius of convergence. 

Since p^ goes to zero, the function ^{A.) must diverge for a finite value of A. Below, we choose A= Ihy convention. 

Remark. In the case of real functions, when the relevant zeros are complex it is often convenient to choose minima 
of the polynomials P^, which satisfy 



choosing, in general, the largest zero for which P^ is small. Other mixed criteria involving a combination of Pjc and 
f ^ can also be used. Indeed, the approximant is not very sensitive to the precise value of p^, within errors. Finally, 
Pk+iiPk) gives an order of magnitude of the error. 

3.2. Functions analytic in a cut-plane: Heuristic convergence analysis 

Although some rigorous convergence results have been obtained [10], there are not optimal. Therefore, we give 
here heuristic but quantitative arguments that show the nature of the convergence of the ODM method. Following 
[20], to simphfy we consider a real function analytic in a cut-plane with a cut along the real negative axis (figure 1) 
and a Cauchy representation of the form 



f'\z) = YjP((PkU'(Pk^z) with Pk(Pk) = 0. 




P'M) = o, 




but the generahzation is simple. Moreover, we assume that 



A(g) (x /e' 



A > 0. 



(10) 



The function E(g) can be expanded in powers of g: 




The assumption (10) then impUes a large order behaviour 



Ek oc i-Ay'Tik -b)~ 



i-Arh-''-^k] , 



exactly of the form displayed in section 1.1. We introduce the mapping 



A 



a > I. 



(11) 



g=P 



ii-Ar' 



The Cauchy representation then can be written as 
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where R{A) is a sum of contributions from cuts at finite distance from the origin. We expand 

E{g(A))^YjPk(p)mt (12) 

k 

with 

Pk(p) - - I d/lA(g(/l))/l * ^ + finite distance contributions. (13) 

For A: -> 00, the factor /t~* favours small values of A but for too small values of A the exponential decay of A{g(A) takes 
over. Thus, Pk(pk) can be evaluated by the steepest descent method. With the Ansatz that at the saddle point A < is 
independent of k and 

Pk ~ R/k, R> 0, 

which implies giA) — > 0, A(g) can be replaced by its asymptotic form (10) for g ^ 0_. At leading order, the saddle 
point equation reduces to 

A(^(l_,r_ln|,|) = 0. (14) 
In what follows, we set R/A = fi, since this is the only parameter. The equation can be rewritten as 

^+l(l_^)-i((^_l)^+l) = 0. 

If the mapping (11) does not cancel all singularities (and this excludes the case of the integral of section 4), then 
PkiPk) cannot decrease exponentially with k. This implies another equation 

-(\-Ar-n\m = o. (15) 

A 

This is indeed the region where the contribution coming from the cut at the origin and from the other finite distance 
singularities are comparable and where the zeros of Pkip) can he. 

Returning to the expansion (12), at g fixed, from the behaviour of pk we infer 

1 - ^ ~ {RIkgfl" ^ A* ~ ^-k'-'"W^)"' . 

In a generic situation, we then expect PkiPk) to behave Uke 

PkiPk) = 0(eC'=""") 

and the domain of convergence depends on the sign of the constant C. For C > 0, the domain of convergence is 

|^|</?C-«[cos(Arg^/a)f. 

For a > 2, this domain extends beyond the first Riemann sheet and requires analyticity of the fiinction E{g) in the 

corresponding domain. 

For C < 0, the domain of convergence is the union of the sector |Argg| < najl and the domain 

\g\>R\Cr[-cos{Mggla)f. 
Again for a > 2, this domain extends beyond the first Riemann sheet. 

3.3. Examples 

For a = 3/2, combining equations (14) and (15), one finds 

p = 4.031233504 , A = -0.2429640300 . 
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Table 2: 



a 


3/2 


2 


5/2 


3 


4 




4.031233504 


4.466846120 


4.895690188 


5.3168634291 


6.1359656420 


-A 


0.2429640300 


0.2136524524 


0.1896450439 


0.1699396648 


0.14003129119 



For a = 2, equation (14) becomes 

A^-fxA-l =0 ^ A=l(jx - + 4) 

For fi = 3.017759126 . . . one recovers the exponential rate of convergence (16). 

In the case of additional singularities, with the additional equation (15), one obtains 

fx = 4.466846120 . . . , A = -0.2136524524 .... 

To give a few other examples, again combining equations (14) and (15) one finds the results displayed in table 2. 

4. Application: The simple integral d = 

For r = 1, the integral (2) in the case d = reduces to the simple integral 

S2nJ 

and the convergence of the ODM method can be studied analytically. 

4.1. The optimal mapping 

Analytic properties suggest that the optimal mapping is given by setting 

Then, / has an expansion of the form 

f{A) = Y,Pk{pmg)f- 

Convergence can be studied analytically. First, / has the representation 



with 

Pk(p)=l:-^ f d^e-^/V/2-p//24)*. 

Setting 

s^/2 = kt, p=R/k, 

one can rewrite the expression as 



Table 3: 

ODMfor the integral d = for g ^ oo: Z(g) ~ g'^'^iMT) * 24^/4 * ^JnlYi^jA). We define [Z{g)g^l'^l 



k 


5 


10 


15 


20 


25 


30 


1/p 


1.131726 


2.35036 


3.34050 


4.5594 


5.5495 


6.8614 


-5 


5.7 * 10-^ 


2.5 * 10-^ 


3.7 * 10-*' 


2.2 * 10-** 


3.4 * 10-'' 


5.1 * 10-" 


ln|<5| 


-5.1578 


-10.5921 


-12.5008 


-17.5923 


-19.4855 


-23.6818 


k 


35 


40 


45 


50 


55 


60 


1/P 


7.7586 


8.9778 


9.9678 


11.1869 


12.1769 


13.3958 


-6 


3.5 * 10-12 


2.5 * lO-i'i 


3.9* 10-15 


2.9 * 10-" 


4.5 * 10-1^ 


3.4 * 10-20 


ln|5| 


-26.3535 


-31.2859 


-33.1625 


-38.0643 


-39.9364 


-44.8208 



For oo, the integral can be evaluated by the steepest descent method. The saddle point equation is 

Rf-ie-{l+RI3)t+l=Q 

and, thus, 

t= - (\+RI3± Vl +/?V9). 
For k odd, the zero corresponds to a cancellation between the two saddle points. This yields the equation 

As expected, one finds 



Pk~T with R = 4.526638689 . . . (R/A = 3.017759126 . . .). 

k 

The minimum is given by one of the saddle points 

Pi(p) oc e-3*/^ = (0.5154353381 . . .f. (16) 
At g fixed, A converges to 1. More precisely, 

A=l- ^jRjki + Oil/k) ^ / ~ e- . 

The approximants converge geometrically on the entire Riemann surface, a situation possible only because the func- 
tion Z{g) has no other singularity at finite distance. 

Numerical verifications. With about 60 terms, the slope is found to be l/fcp^ ~ 0.2209 in agreement with the 
prediction \IR - 0.2209 (once even-odd order oscillations are taken into account). The logarithm of the error has a 
slope 0.696/0.685 to be compared with the prediction 3//? = 0.66 (see table 3). 

4.2. An alternative mapping 

Another mapping that also regularizes the point at infinity is 

A 

g=P 



a -Ay 
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Table 4: 

ODMfor the integral d = Ofor g = 5: we define [Z(g)]exact - [2(g)]oDM = S 



k 


5 


10 


15 


20 


25 


30 


1/p 


0.5918 


1.0297 


1.5627 


2.0779 


2.5865 


3.1376 


1^1 


1.1 * 10-^ 


3.7 * 10-5 


1.7* 10-*" 


9.2 * 10-^ 


1.1 * 10-^ 


3.5 * lO-"* 


ln|<J| 


-6.7454 


-10.2069 


-13.2837 


-13.8898 


-16.0103 


-19.4614 


k 


35 


40 


45 


50 


55 


60 


1/P 


3.6167 


4.1877 


4.6557 


5.3021 


5.6959 


6.2458 


-5 


3.4 * 10-'' 


8.0 * 10-"' 


1.0* 10-'" 


2.9* 10-" 


1.2* 10-" 


2.4 * 10-'2 


\n\6\ 


-19.4706 


-20.9453 


-22.9796 


-24.2450 


-25.0907 


-26.7433 



Numerical results for ^ = 5 are displayed in table 4. Finally, for k = 65, 70, 

l/Pk = 6.7479 , 7.5724 , In |5| = -30.3657 , -28.9505 
On the average between k = 5 and k = 70, 

Pfe~/?/A; with/? = 9.75. 
The results displayed in table 2 lead to the prediction 

R = 9.2039 . 

The error is about 

p/4 

5 = -1.59-r-r 

and 1.59 has to be compared with the expected asymptotic value 1.74 if one assumes convergence for all g > 0. 

5. The quartic anharmonic oscillator: d = 1 

For r - I, the path integral corresponds to the quantum Hamiltonian 

rj _ 1 „2 , 1 ^ , ^ ^4 

H - +^xf + —x. 

The eigenvalues E of H are given by the solution of the time-independent Schrdinger equation 



where tlf(x) is a square -integrable function. 

As an example, we consider the perturbative expansion of the lowest eigenvalue, the ground state energy. Varia- 
tional arguments and scaling suggest the mapping 

Then, 

k 

9 



Large order behaviour (section 1.1) and a steepest descent evaluation (table 2) lead to the prediction 

Pk ~ R/k with R = nA = 32.25 .... 

Then, A converges to 1 as 

a = 1 _ / + o{k-'l') ^ A* ~ e-«'"*"'/«''' with /?2/3 = 10.131 .... 

An unbiased fit of the numerical data for A: < 60 yields results within 10% of the predicted values. Finally, a fit of the 
relative error for g — > oo yields [20] 

P..i(p.)oce-'-^*"'. 

The relative error at order k is thus of order e~*"'('-^+*^''« One finds convergence for 

0.95 + |^r^/^ cos(|Arg^) >0. 

The corresponding domain contains a section of the first Riemann sheet and extends to the second Riemann sheet for 
1^1 large enough. 

6. field theory iad = 3 dimensions 

In [20], the ODM method has been applied on functions of the initial parameter g of the action (1) rather the renor- 
malized parameter introduced in section 2. Then the point of physical interest is ^ — > oo, which corresponds to the 
zero g* of the ;6-function (8). Due to UV divergences, a needed regularization and renormalization, scaling arguments 
are no longer applicable to determine an appropriate mapping. The relation (9) between initial and renormalized 
parameter shows that, for ^ ^ oo, physical observables have an expansion in powers of g~", where the exponent 
(X) - /3'Q*). This then suggests the mapping 

A 



8=P 



(1 



(18) 



but the difficulty is that a) has to be inferred from the series (8) itself. The results obtained in this way [20] are 
consistent with those obtained in [12] {co = 0.80(1) from Borel transformation and mapping), but empirical errors are 
more difficult to assess. Also the expected rate of convergence is of order e^'""^* *' " = g-const.*"^^ which is rather slow 
(see table 5). See reference [20] for details. Finally, the information about the large order behaviour cannot easily be 
incorporated. 



Table 5: 



Series for the exponent co summed by ODM in the (pi field theory with tJi^ — 0.79 and dwcai./dwin — -0.6. 



k 


2 


3 


4 


5 


6 




0.552 


0.754 


0.711 


0.767 


0.759 



Here, to illustrate the flexibility of the method, we work directly with functions of g. We also take into account 
the covariance of the jS-fimctions under a change of parametrization: 



(19) 



This transformation law is such that the derivative of the yS-function at a zero (a fixed point), which is a physical 
observable, remains unchanged. 
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Equation (19) suggests a mapping of the form 



with a suitable choice of the parameter a, since unlike a mapping of the form (11), it introduces no new singularity. 
We thus set 

M^) = - — —Km). 

ap 

A few trials, without trying to optimize, suggest the value o; = 3/2 and this is the value we have adopted. The results 
for the zero g*, for comparison with the method outlined in section 2, and the exponent u) - /3\(A*) are given in table 
6. The order of magnitude of the errors can only be estimated by the sensitivity to the precise choice of the parameter 
p (zero of last term or its derivative, for example). The indications are Ag* » Aa> » 0.006. In the case of complex 
zeros, we have given only the real part in the tables. 



Table 6: 

Results for the zero g* of the ^-function and the exponent oj summed by ODM in the (pi field theory. With the method 
of section 2 ([12], one finds |* = 1.411 ± 0.004 and oj = 0.799 ± 0.011. 



k 


3 


4 


5 


6 


7 


r 


1.09871 


1.39330 


1.41771 


1.41737 


1.41744 


0) 


1 


0.7984 


0.7804 


0.7806 


0.7807 



Other critical exponents are obtained from the two RG functions 

y-i(|) = l-i|+i|2-0.0230696212|3+0.0198868202|4- 0.0224595215 15 
+ 0.0303679053 f - 0.046877951 g'^ + 0{g^), 
?7(|) = 0.0109739368 f + 0.0009142222 f + 0.0017962228^'^ - 0.0006537035 f 
+ 0.00127491001^ - 0.001697694 1'^ + 0{g\ 

by setting g - g* ■ We have also summed the series for y~'^{g) and ri(g)/g^, and independently the series for the 
function v~^(g) although it is related to the y and 77 by y(g) = v(g)(2 - T](g)). A verification of this relation after 
summation gives an indication about the errors. The results, displayed in table 7, can be compared with the results of 
[12]: 

y= 1.2396 ±0.0013, v = 0.6304 ± 0.0013 , 77 = 0.0335 ± 0.0025 . 
One notices a reasonable consistency. 



Table 7: 

Critical exponents y — y{g*), v — v(g*) and rj — i](g*) for g* — 1.411. 



k 


3 


4 


5 


6 


7 


7 


1.23717 


1.23486 


1.23845 


1.23820 


1.23923 


V 


0.62521 


0.62486 


0.62746 


0.62771 


0.62865 


J] 




0.0290 


0.0289 


0.0297 


0.0306 
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7. Convergence proofs and physics applications 

Since the ODM summation method has been proposed, convergence has been proved in specific cases [10], how- 
ever, to our knowledge, systematic mathematical investigations are still lacking and would be most welcome. 

Notice that in [8J a special case of the ODM method has also been proposed. Moreover, methods known as linear 
or scaled delta expansion, or optimized perturbation theory, introduced later, often reduce to special cases of the ODM 
method (see, for example, [9], [4]). 

The ODM method has also found a number of useful apphcations in physics (sometimes under difi'erent names). 
In [18], the problem of the hydrogen atom in strong magnetic fields has been considered. The summation of the weak 
field series expansion has led to precise determinations of the ground state energy for very strong fields. In [13], 
the location of the Bender-Wu singularities of the quartic anharmonic oscillator has been determined numerically 
(for the anharmonic oscillator see also [2]). In [11], the equation of state of physical systems belonging to the Ising 
model universality class has been determined numerically by a combination of Borel transformation and mapping and 
the ODM summation method. Another application has been the Bose-Einstein condensation problem [21], [14] (for 
similar problems see also [6]). Finally, let us also mention the appUcation to the Gross-Neveu model [15], [16]. 
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